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ABSTRACT 

The random-bond Ising model on the square lattice has several disordered 
critical points, depending on the probability distribution of the bonds. There 
are a finite-temperature multicritical point, called Nishimori point, and a zero- 
temperature fixed point, for both a binary distribution where the coupling 
constants take the values ±J and a Gaussian disorder distribution. Inclusion 
of dilution in the ± J distribution ( J = for some bonds) gives rise to another 
zero-temperature fixed point which can be identified with percolation in the 
non-frustrated case (J > 0). We study these fixed points using numerical 
(transfer matrix) methods. We determine the location, critical exponents, and 
central charge of the different fixed points and study the spin-spin correlation 
functions. Our main findings are the following: (1) We confirm that the 
Nishimori point is universal with respect to the type of disorder, i.e. we 
obtain the same central charge and critical exponents for the ± J and Gaussian 
distributions of disorder. (2) The Nishimori point, the zero-temperature fixed 
point for the ±J and Gaussian distributions of disorder, and the percolation 
point in the diluted case all belong to mutually distinct universality classes. 
(3) The paramagnetic phase is re-entrant below the Nishimori point, i.e. the 
zero-temperature fixed points are not located exactly below the Nishimori 
point, neither for the ±J distribution, nor for the Gaussian distribution. 
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1 Introduction 



The problem of disordered magnetic systems has attracted great interest in the past years 
and many questions still remain unanswered. An interesting problem is the universality 
class of second order phase transitions in two-dimensional systems. The random-bond 
Ising model (RBIM) is one of the simplest and best known of these systems [1-4], but 
exhibits rich enough behaviour to give a general understanding of the problem. Further 
interest in the RBIM stems from analogies with the quantum Hall transition [5-8] and 
applications in coding theory [9-12]. 

For a small amount of randomness, the universality class of the RBIM remains un- 
changed, presenting only logarithmic corrections in some correlation functions [1]. For 
some particular kind of randomness, namely dilution, one can show that another non- 
trivial fixed point corresponds to a percolation universality class at zero temperature, 
since it becomes a purely geometric problem of having a thermodynamic number of spins 
within the same cluster. It is important to notice that in this case, all the randomly 
distributed bonds are non-negative. 

The situation is quite different if some negative bonds are allowed in the probability 
distribution. For certain distributions with negative bonds, Nishimori [2, 3] has shown 
that some exact statements can be made about physical quantities. There is in particular 
a line in which the internal energy can be calculated exactly, known as the Nishimori 
line. The interest of this line goes further, since it has been shown that this line is in- 
variant under renormalization group (RG) transformations [4]. Since this lines crosses 
the ferromagnetic to paramagnetic transition line, the intersection point, known as the 
Nishimori point is a fixed point. Examples of probability distributions satisfying the 
Nishimori condition are the Gaussian, ± J binary and 0, ± J with appropriately chosen 
weights (see below). Despite recent analytical approaches [6], the exact characterization 
of the universality class of this non-perturbative fixed point is still unknown. The under- 
lying field theory describing this point is certainly a good representative of the disordered 
fixed-point behaviour with very interesting and rich phenomenology [13,14]. In a previ- 
ous letter [15], we have studied numerically the critical exponents and central charge of 
the Nishimori point in the ±J RBIM. We found in particular that its universality class 
does not correspond to the one of percolation, as one could have imagined in view of 
earlier numerical investigations of this model [16-18]. This is derived, first, from numeri- 
cal estimates for the critical exponents and central charge which differ significantly from 
those of percolation, and confirmed by an analysis of higher moments of the correlation 
functions (which are all equal in percolation and only in pairs for the Nishimori point, 
see [15] for details and later in this paper). 

The conclusion of a different universality class was confirmed in another numerical 
analysis by Merz and Chalker [19]. Thanks to a mapping to a network model, these au- 
thors reached big lattice sizes and high accuracy in the measurement of critical exponents. 
An interesting remark made by these authors concerns the dual theory of the RBIM 1 , 
in which the different moments of the disorder field acquire negative dimensions [20]. 
More recently, Nishimori and Nemoto [21] used a generalized model with self-duality, 
and conjectured that the projection onto the RBIM gives the phase boundary for this 

1 Because of the randomness introduced in the bonds, the model is not self-dual as the pure Ising 
model. 
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Pc 


Method 


Reference 


Nishimori point in the ±J model 


0.111 ±0.002 


Transfer matrix 


[30] 


0.114 ±0.003 


Series expansion 


[17] 


0.1128 ±0.0008 


Non-equilibrium 


[31] 


0.1095 ±0.0005 


Transfer matrix 


[18] 


0.1094 ±0.0002 


Transfer matrix 


[15] 


0.1093 ±0.0002 


Fermionic transfer matrix 


[19] 


0.110028 


Duality 


[21] 


< 0.178203 


Rigorous upper bound 


[32] 


T = critical point in the ±J model 


~ 0.099 


Series expansion 


[22] 


0.105 ±0.01 


Matching algorithm 


[23] 


0.095 <p c < 0.108 


Matching algorithm 


[24] 


0.104 ±0.001 
0.106 ±0.002 
0.115 


Exact ground states 
Ground state enumeration 


[25] 
[27] 


0.1031 ±0.0001 


Exact ground states 


[12] 


0.103 ±0.001 


Exact ground states 


[28] 



Table 1: Overview of estimates for p c at fixed points in the two-dimensional ± J random- 
bond Ising model. The first part of this table is for the Nishimori point, the second part 
for the zero-temperature critical point. 

model. This result permits in particular to locate the Nishimori point analytically. The 
conjectured location is however outside the accuracy range of the most recent numeri- 
cal works on the ±J binary disorder case [15, 19] and the validity of this conjecture is 
certainly a very interesting open issue. 

The pure and the Nishimori point, when present in the phase diagram, are not the only 
non-trivial fixed points of the model: there is in any case a zero-temperature fixed point. 
While Nishimori's results state rigorously that this point cannot be located at a higher 
concentration of "impure" bonds, analytical and numerical works on the ±J model [12, 
22-28] tend to conclude that it is located at a smaller density of impurities, indicating a 
re-entrance of the ferromagnetic phase. The properties of the zero-temperature point vary 
considerably with the kind of disorder introduced. For a symmetric Gaussian distribution, 
it has been shown [29] that the lowest energy configuration is unique (modulo the Z2 
symmetry) with probability one. One is tempted to check the extension of this result to 
our zero-temperature point, and check for example that, for a given configuration of the 
disorder, any spin-spin correlation function is 1 or —1. Then, all the odd moments of the 
spin-spin correlation functions are equal, and all the even moments are just equal to 1, 
a result that is similar, but not exactly identical to the percolation case. The situation 
is much more subtle for distributions like the ±J one, since frustration plays a crucial 
role. For a generic configuration of disorder, the lowest energy states are expected to be 
highly degenerate. The results obtained by different techniques for the location of both, 
the Nishimori point and the zero-temperature critical point are summarized in Table ^ 
for the ±J model. As can be seen from the most recent results, the zero-temperature 
critical point seems to be located at a concentration of "impure" bonds strictly smaller 
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Figure 1: Phase diagram of the two-dimensional ±J random-bond Ising model. The 
arrows represent the flow under the renormalization group. 

than the one of the Nishimori point. A schematic picture of the phase diagram of the 
± J case in the p — T plane (p being the number of antiferromagnetic bonds) is shown 
in Fig. ^ where the separation of the paramagnetic and ferromagnetic phases, as well 
as the Nishimori line are drawn. The location of the Nishimori point N and the zero- 
temperature critical point are representative of the results shown in Tabled For the case 
of the Gaussian disorder, a similar diagram may be obtained by replacing the parameter 
p with the variance of the distribution of disorder, a, although the shape of the Nishimori 
line is different. 

Given a simple phase diagram as shown in Fig. ^ it is tempting to trace the massless 
flow backwards from the pure Ising model at p = to the Nishimori point N [33,34]. 
However, we will argue in this paper that in a generalized parameter space there are in 
fact several fixed points on the massless surface with different universality classes such 
that a backtracking procedure is in general ambiguous. Indeed, the problems encountered 
in [34] were one of the original motivations for our earlier numerical study of the Nishimori 
point [15] which we continue and expand here. 

In this paper we provide an extensive numerical analysis of the random-bond Ising 
model for three different kinds of probability distributions. We first consider in parallel 



3 



the case of Gaussian and ± J distributions for which we estimate the location of the fixed 
point along the Nishimori line using domain-wall techniques. We compute the free energy 
with a transfer matrix technique [35] and obtain the central charge by analyzing the finite- 
size corrections. We also compute the spin-spin correlation functions on long strips to 
obtain the magnetic exponent. This and other measurement to check consistency clearly 
show that the models with ±J and Gaussian distributions share the same universality 
class at their respective Nishimori points. The compatibility of the results between the 
± J and Gaussian cases further elucidates the nature of the universality class of this point 
which was recently argued to be different from the one of percolation [15, 19]. We next 
consider dilution, which can be modeled by a probability distribution allowing the values 
of ± J and for the coupling constants. In the case of pure dilution (allowing only values 
+J and for the coupling constants), the other fixed point apart the one of the pure 
model is located at zero temperature and dilution q c = 1/2 and simply corresponds to 
bond percolation. The ± J, distribution allows us to study the behaviour in the critical 
line between the Nishimori point and the percolation point. In particular, by studying 
the value of the effective central charge for different strip widths and its evolution towards 
larger sizes we show that the percolation fixed point is repulsive along this critical line in 
favour of the Nishimori point. The same technique is used to confirm that the Nishimori 
point is unstable with respect to the pure Ising fixed point when moving on the critical line 
connecting these two points. We finally address the problem of the zero-temperature fixed 
point. Although for the case of the ± J distribution there is an extensive list of numerical 
works indicating re-entrance of the ferromagnetic phase (see Table ^), to our knowledge 
there was no conclusive evidence of the same fact for the Gaussian distribution. We show 
here clear evidence for the re-entrance of the phase also for the Gaussian distribution. 
We compute also different moments of spin-spin correlation functions and the magnetic 
exponent. Our results in the magnetic sector clearly show that the universality class 
of this zero-temperature critical point is once again different from the one of the finite- 
temperature Nishimori point as well as percolation, in contrast to what one could have 
thought considering previous numerical results in two dimensions [17,25]. 

2 Some definitions 

In this section, we present some definitions which will be used throughout this paper. 
We will consider two kinds of probability distributions P{J) for the bonds: a discrete 
distribution where the coupling constants can take values ±1 and 0, and a continuous 
Gaussian distribution. In the case of a discrete distribution, two subclasses can be 
considered, the first one allows only the values ±1 while the second and more general 
one, allows also the value for the coupling constants. This last case corresponds to 
dilution. In any case we can imagine a phase diagram in which the vertical axis is given 
by temperature and the horizontal axis by a parameter representing the strength of the 
disorder (see Figure [TJ). 

Let us define the Hamiltonian of the two-dimensional Ising model to be: 



were S(i) = ±1 and (i, j) means nearest neighbours on a square lattice. The variables 
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Jij are random and independently chosen with a probability distribution function P(J). 
As usual, the Kronecker delta function is zero if the spins are different and one if they 
are equal. To establish the relation to other conventions for the energy note that it can 
be expressed in terms of products of Ising spins via Ss(i),s(j) = S(j) + 1) /2. 

The Nishimori line is defined by the condition [2]: 

P{—J) — exp(-/3 J) P(J) , (2.2) 

where (3 = -^j, is the inverse temperature (from now on we choose the convention ks = 1)- 
Let us analyze the Nishimori condition for the different distributions considered here 

• J = ±1: The distribution is characterized by the concentration of antiferromagnetic 
bonds p: 

P{J) = {l-p)S{J-l)+pd{J + l) . (2.3) 
Eq. (|2.2|) gives the following condition for the Nishimori line 

= In , (2.4) 

which is schematically depicted in Fig. ^ The line extends from the (attractive) 
fixed points given by T = 0, p = and T = oo, p = 1/2 and crosses the Para- 
Ferro transition line at the critical concentration p c . In the first part of Tabled we 
summarize estimates for the location of the critical point p c . With the exception of 
a conjectured duality property [21] and a rigorous upper bound [32], all the other 
results come from numerical simulations. 

• Gaussian distribution: 

™ = S ("^) • (2 ' 5 > 

with the following condition 

for the Nishimori line. In the following, we will choose Jq — 1 without any loss of 
generality. Thus, the distribution is characterized by the parameter a along the 
Nishimori line which extends from the (attractive) fixed points given by T = 0, 
<T = and T = oo, a = oo and crosses the Para-Ferro transition line at the value 
a c . The first numerical characterization of the critical point was given in [16], with 
a value of a c ~ 0.97. 

• Finally, the binary distribution can be generalized to include dilution. In this case, 
we write: 

P(J) =q5( J) + (l-q-p) S( J- l)+pS( J +1) . (2.7) 

The case q = corresponds to the binary distribution discussed above, while p = 
gives the ferromagnetic diluted model. In the diluted model there is only a 
fixed point located at zero temperature and q c = 1/2, apart from the pure model 
fixed point [36]. This T = fixed point corresponds to percolation. The reason 
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Figure 2: Phase diagram of the diluted ±J, random-bond Ising model, p and q are 
the concentration of negative and zero bonds, respectively. The arrows represent the flow 
under the renormalization group. The dashed line represents the intersection between 
the Nishimori surface and the Ferro-Para transition surface. 
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for this is that in the absence of frustration all the spins belonging to the same 
cluster must point to the same direction at zero temperature. Then, whenever 
there is percolation, or a cluster containing a macroscopic number of spins, there 
is a macroscopic magnetization, and this happens until the critical concentration 
of bonds 1/2 (see for example [36,37]). 

The Nishimori surface in the T-p-q space is now given by: 

^ln(i^P). (2.8, 

The intersection of this surface with the ferromagnetic-paramagnetic transition 
surface gives a renormalization group invariant line [4] whose end points are the 
percolation q = 1/2, T = fixed point and the finite-temperature Nishimori point 
N of the purely binary depicted in Fig. [TJ In the T — p — q phase diagram, 

the location of the Nishimori point, here called N', must be within the intersection 
line mentioned above (the dashed line in Fig. EJ). In this sense, the point N depicted 
in Fig. Q is the representative in the T — p space of the more general location of the 
Nishimori point denoted by N' in Fig. El 



3 Domain-wall free energy 

In the following section we will discuss the domain-wall free energy in a manner very 
similar to [16]. For a strip of with L the domain- wall free energy dz is defined as 

d L = I? (f? - f { L a) ) , (3.1) 

where is the free energy per site = ln L z ^f of a strip of width L and length M with 

periodic boundary conditions 2 and = the corresponding one with antiperiodic 

boundary conditions. 

g?£ measures the free energy associated to a domain wall in the system. We will first 
consider the ±J distribution of disorder. For fixed parameters p (/3) in the disordered 
(paramagnetic) phase, one should have lim.L-^ood>L — > while dL should diverge in the 
ordered phase (d^ — > oo as L — > oo). At the fixed point p c (/3 C ), g?l should converge 
quickly with L. We can therefore use crossing points between d^ and d^ 2 to obtain 
finite-size estimates for p c . 

Let us consider first the ± J distribution of disorder, and f L °^ are computed with 
the transfer matrix technique (see, e.g., [35]) on strips of length 10 6 . Averages over up to 
iV m 4000 samples of such L x 10 6 strips are taken in order to average over randomness 
and to determine the statistical error. It is useful to fix the number of bonds on each 
L x 10 6 strip to approximate the chosen value of p as accurately as possible: In our 
implementation such a strip has L ■ (2 • 10 6 — 1) bonds out of which we select the integer 
closest to pL ■ (2 • 10 6 — 1). We found that this method leads to 5 times smaller error bars 
for f L p ^ than if one selects each bond separately at random with probability p, i. e. for the 
same precision one needs 25 times less samples with this method compared to generating 

2 Note that our sign conventions differ from the standard ones. 
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each bond at random with probability p without constraining their total number. For 
di the improvement is even bigger and the error bars become 10 times smaller - or one 
needs 100 times less samples for the the same accuracy. 
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Figure 3: Position p of crossings between d^{p) and di 2 (p) on the Nishimori line of the 
±J random-bond Ising model with L 2 < 12. 

Next we will present some details of the analysis leading to the results presented in [15] 
- a very similar computation was also performed in [19]. Fig.OJshows a doubly logarithmic 
plot of the crossing point p c (Li,L 2 ) between d^ip) and dL 2 (p), using p c = 0.1094 in 
the shift of the vertical axis. It is known for the pure Ising model that the finite-size 
corrections to p c scale roughly ~ L~ 2 with L = (Li + L 2 )/2 [38]. In order to extrapolate 
p c to the thermodynamic limit, one can therefore first use a plot of finite-size estimates for 
p c as a function of L~ 2 and extrapolate to the vertical axis. A reasonable extrapolation 
with a generous error bar is 

p c = 0.1093 ±0.0004. (3.2) 

Even if the correction is probably not of the form 1/L 2 in the present case, one can see 
in Fig. El that the finite-size correction to p c is well described by a power in L. For an 
improved extrapolation we therefore use the following form for a fit: 

p c (L 1 ,L 2 ) = Vc + ot\ Ll \ L ) ■ (3.3) 



We then find [15] 

p c = 0.1094 ±0.0002 (3.4) 
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and an exponent 



f = 1.5 ±0.3. 



(3.5) 



The two estimates ()3.2j) and ()3.4j) agree well with each other - the error bar of the second 
one is just a little smaller. The exponent ()3.5|) cannot be determined very accurately. 




Figure 4: Slope A(L) of the domain-wall free energy for the ±J random-bond Ising 
model as a function of 1/L. 

The estimates (|3.2j) and (|3.4j) should be compared with other estimates which are 
summarized in Tabled At this point comparison with the zero-temperature transition is 
meaningful only if one assumes that the T = transition point is located exactly below 
the Nishimori point. Our estimate agrees within error bars with other transfer matrix 
computations [18,19,30]. However, our estimate falls outside error bars of estimates 
obtained by other methods: It is smaller than the estimates of [17, 31] (in the latter 
case significantly). This discrepancy becomes clearest if one notices that several of our 
finite-size crossings have already reached the region p pa 0.1097 (see also Fig. 2 of [15]). 
If one now assumes that these crossings decrease monotonically with increasing L (as the 
data in Fig. El indeed does), one does not need to actually carry out the extrapolation 
and estimate its error in order to conclude that our extrapolated value for p c must end up 
below the error windows of [17,31]. At this point it is particularly reassuring to observe 
that our result is fully consistent with the one of [19] which has been obtained using also 
substantially wider strips (up to L = 64). Furthermore, (|3.4j) is in complete agreement 
with other quantities to be discussed in later sections. Thus, we observe a significant 
difference between the analytic estimate of [21] obtained from a duality argument and 
our value for p c as well as the one of Merz and Chalker [19]. 
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Now let us look at a first critical exponent, namely the correlation length exponent v 
along the Nishimori line. If one assumes the scaling form 

d L (p- Pc ) = d((p-p c )L^) , (3.6) 

one can expand around p c and finds 

dh{p — Pc) ~ const. — A(L) p (3.7) 

with 

A{L)~I}/ V . (3.8) 

Then one can fit d>L close to p c by a linear function and extract A(L). 

Fig. |3] shows the values for A(L) determined in this manner on a doubly logarithmic 
scale 3 . One can see that they follow indeed a power law. Using ()3.8j) we extract 

v= 1.48 ±0.03, (3.9) 

which amounts to a slight correction of the value v = 1.33 ± 0.03 obtained in [15]. The 
result ()3.9j) is now in excellent agreement with the value v = 1.50 ± 0.03 obtained by 
a fermionic transfer matrix [19]. The most recent results v « 1.5 do no longer agree 
well with v = 1.32 ± 0.08 obtained by high-temperature series [17]. However, we already 
observed above that the series expansion method does not yield a very accurate estimate 
for p c either. Furthermore, the value v = 4/3 characteristic for percolation (see e.g. [37]) 
now falls outside numerical errors and thus it is possible to conclude that the Nishimori 
point is not in the universality class of percolation already on the basis of the exponent 
v. Note that the result (|3.9|) does not involve locating p c precisely and should therefore 
be independent of errors which may have been made in the location of p c . 

Figures with scaling collapses of d^ were presented in [15,39] with v = 1.33 and in [19] 
with v = 1.50. While they verify that di obeys indeed the scaling form ()3.6|) . the fact 
that reasonable collapses can be obtained for different values shows that such a collapse 
is not a good criterion for determining v. 

We now present the same domain-wall free energy analysis for the Gaussian distribu- 
tion of disorder. Again, and are computed with the transfer matrix technique 
for small sizes. For larger sizes (L = 15,18,21,24), we employed a different algorithm, 
developed by Merz and Chalker, which uses a mapping to a network model, see [19] for 
details. For each point we averaged over 10 000 samples of size L x 10 6 . In Fig. we 
show a plot of the crossing point a(Li,L 2 ) between d^^cr) and o?l 2 (<t). Contrary to the 
±J case where the crossing was converging quickly, here we see that the situation is 
much more complex. For the smaller sizes with L 2 = L\ + 2 and 2/(Li + L 2 ) > 0.15, it 
apparently converges to o ~ 0.9815. For larger sizes, we see a different behaviour. We 
explain this change of behaviour by the existence of an irrelevant operator. If such an 
operator exists, we expect that close to the critical point, the domain- wall free energy 
can be expanded as [16,19,38]: 

d L {a) =a + b{a -a c )L 1/u + cL~ x . (3.10) 

3 Our data for di, is in perfect agreement with that of [19] where we overlap. However, the windows 
in p used for the present estimates of A(L) differ from those used in [15] leading to slightly different 
results. 
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Figure 5: Position o of crossings between d^a) and di 2 (a) on the Nishimori line of the 
Gaussian random-bond Ising model. The figure also contains a linear fit to the large size 
data and the value a = 1/1.02177 predicted by duality arguments [21]. 
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Figure 6: Slope of the domain-wall free energy of the Gaussian random-bond Ising model 
function of 1/L. 
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Figure 7: Domain-wall free energy of the Gaussian random-bond Ising model at a c = 
0.97945 as a function of L. We also show a plot of the fit to the form (j3.10|) with L > 7. 



Here we have neglected a further contribution coming from the thermal exponent 
(see [19]) which is sub-leading compared to the (a — a c ) term. The last contribution 
corresponds to an irrelevant operator, of dimension 2 + x, associated to the Nishimori 
point. In principle it should also be present in the scaling of the Nishimori point with 
the ± J distribution although much bigger sizes may be necessary for its observation. If 
we take only the points with 2/(L\ + L 2 ) < 0.1 then the data in Fig. EJseems to converge 
linearly to a point that we determine by a fit to be 

a c = 0.97945 ± 0.00004 . (3.11) 

In Fig. El we also show the plot corresponding to this linear fit (for the larger sizes) as 
well as the value a = 1/1.02177 = 0.978694 predicted by duality arguments [21]. While 
the two values are close, they are still not compatible from our analysis. If one changes 
the parameters of the linear fit, allowing for smaller size data to be taken in account, the 
smallest possible value that we get is a = 0.97939 ± 0.00002 which is still not compatible 
with the duality argument value. The only other numerical estimate of a c that we are 
aware of is the one of McMillan [16], a c ~ 0.93 — 0.97. This result is not expected to be 
highly accurate since it is obtained on small systems of linear sizes up to L = 8. 

Next we turn to the determination of v. Again, we compute it by looking at the 
slope of the domain- wall free energy close to the critical point. In Fig. |Ul we show the 
slope as a function of the size. The slope is determined in the following way: for each 
size, we computed d^a) for a = 0.980 and o = 0.982. (For the sizes that we considered, 
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the effective critical value of cr c (L), determined as the crossing of the domain-wall free 
energies, is always in this range, see Fig. EJ). Thus the slope is determined directly from 
the difference between these two values for each size. We also show in this figure the best 
fit to the form ~ L 1 ^. This fits perfectly the data for L > 4. We obtain v = 1.52(3) 
which is very close to the result for the ±J disorder case. Note that using a transfer 
matrix method McMillan [16] obtained v « 1.58 for the Gaussian random-bond Ising 
model. 

In Fig. we plot di,{a) for the critical value a c = 0.97945 determined previously. 
One expects then that the non-constant part comes entirely from the leading irrelevant 
operator. A fit to the form (|3.10|) . also shown in Fig. [7[ gives excellent results if we 
remove the data with L < 6. We obtain an exponent x = 0.43(3). Thus we predict the 
existence of an irrelevant operator of dimension 2.43(3). We observed previously that 
a c (L) converges linearly (in l/L) towards a c {L — ► +oo), see Fig. EJ. A simple calculation 
starting from eq. (|3.10j) shows that the correction is of order L~^^ u+X ^ and the linear 
correction corresponds to x ~ 0.33. This is rather close to the numerical result, the 
difference being easily explained by taking into account further irrelevant operators. 



4 Free energy and central charge 

In this section, we use the free energy to deduce the central charge. We obtain results 
for both the ±J and the Gaussian distribution of disorder and compare these results to 
check universality. 

As a byproduct of the measurements of the domain-wall free energy presented in 
the previous section, one also obtains the free energy . We will first discuss the ±J 
disorder case. Averaging has been performed over sufficiently many samples to obtain 
statistical errors of typically 5f^ < 4- 10~ 6 . Note that due to the aforementioned round- 
off error of one bond, p has an error of about 10~ 7 on the strip sizes considered which is 
not substantially below 5f^\ 

The (effective) central charge c can be estimated 4 from this data since it appears as 
the universal coefficient of the first finite-size correction [40] 

& = & + £ + ("J 

Note that the leading term fso is not universal. The universal l/L 2 term originates from 
the energy- momentum tensor T which has scaling dimension 2. The next correction is 
expected to arise from T 2 and should therefore give rise to a term of the form L~ 4 

& = & + £ + £ + (4.2) 

Estimates for c are obtained by fitting the free energy data in intervals L < L < L max to 
()4.1|) or (|4.2|) . We have considered only intervals with L Q < L max — 3 (at least four-point 

4 For a non-unitary theory what we call 'c' is in fact c c ff = c — 12 5 m i n where c is the central term 
in the operator product expansion of T with itself, and <5 m i n the smallest dimension contributing to the 
free energy. See section UTTl for further details. 
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fits). For a given form of the fit and interval of lattice sizes, we have chosen c to lie in 
the center of the resulting fits and adjusted the error to include all fits. The estimates 
for c obtained in this manner are shown in Fig. |SJ Fits to (|4.1j) for lattice sizes 5 < L < 9 
are denoted by diagonal crosses and those for 6 < L < 12 by plusses. Fits to ()4.2|) were 
performed for 4 < L < 6 with L max kept fixed at the largest available system size. The 
corresponding estimates are denoted by boxes. 
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Figure 8: Estimates for the central charge c along the ±J Nishimori line (|2.4jl in the 
vicinity of the Nishimori point obtained from (|4.1|) for 5 < L < 9 ('x'), 6 < L < 12 ('+') 
or (|4.2|) with 4 < Lq < 6 (boxes). The grey shaded areas denote our final confidence 
intervals for p c (vertical) and the central charge (horizontal). 

One notices that the estimates in Fig. |H] are almost independent of p in the range 
considered. On the one hand, this means that the estimates for c are stable. On the 
other hand, this implies that estimates for the central charge are not a good tool for 
locating the critical point precisely for the ±J distribution of disorder (the situation 
will be different for the Gaussian distribution, see below). Furthermore, one observes 
a slight trend of the fits ()4.1j) obtained with only a L~ 2 term to shift to smaller values 
as the range of system sizes considered is increased while the fits (14. 2j) with a L~ 4 term 
included converge rapidly with system size. Therefore, we quote as a final estimate with 
a generous error bar [15] 

c = 0.464 ±0.004, (4.3) 

which is shown by the grey shaded horizontal bar (the estimate (|3.4|) for the location of 
the critical point p c is shown by the grey shaded vertical bar). In any case, all estimates 
obtained from ()4.1|) with L max = 12 and L > 6 satisfy c < 0.469. Furthermore, also 
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Figure 9: Central charge vs. a along the Nishimori line for the Gaussian distribution of 
disorder, with L max = 8. 

consideration of other possible finite-size corrections to (j4.1j) yields results which are 
consistent with (|4.3j) . Thus, we can safely exclude the value c = 5v ^| n2 ~ 0.47769 for 
percolation in the Ising model [41] even if the absolute difference is not big. 

We now turn to the Gaussian distribution of disorder. In this case, we will argue 
that the central charge, along the Nishimori line, gives a practical way of locating the 
multicritical point as the point where the estimates of the central charge have a maximum. 
We explained above that, for the ± J distribution, the estimate of the central charge was 
nearly independent of p making it difficult to locate the maximum. For the Gaussian 
distribution, the parameter p is replaced with the variance of the distribution a, see 
eq. (J2.5j) . Thus, the situation is different since it is possible to perform the simulations 
with the same disorder configuration for each variance a. More precisely, one defines 
a bond J in the following way: J = Jq + aj, j being the random part obtained from 
a random number generator which we choose to be the same for different a. As a 
consequence, most of the measurements will be correlated and in particular this will be 
true for the free energy and the deduced central charge (as was also the case for the 
domain- wall free energy measured in the previous section). 

The central charge is shown in Fig. El where we present the effective central charge 
obtained by a fit to eq. ()4.1|) using free energies in the range [L Q , L max = 8] for increasing 
L . For the Gaussian distribution, we cannot suppress disorder fluctuations by fixing 
the number of bonds with a certain value globally. As a consequence, the data shown 
in Fig. El which has been extracted from runs of the same size as for the ± J disorder 
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Figure 11: Central charge at the multicritical fixed point on the Nishimori line for the ± J 
and Gaussian distributions of disorder. We also show for comparison finite-size estimates 
of the central charge for percolation. The horizontal line shows the infinite system limit 
c = 5 2 for percolation on Ising clusters [41]. 

(10 000 samples of strips of sizes L x 10 6 ) has much bigger error bars. Still these larger 
error bars are partially compensated by the fact that the measurements are strongly 
sample-to-sample correlated, and thus the maximum is easy to identify: we have a stable 
maximum at a ~ 0.98. In Fig. EH we show the measured central charge close to the 
maximum with L max = 10. An estimate of the maximum gives for all the ranges of L 
employed a value of a c ~ 0.9805 ± 0.0005, with no measurable change of the central 
charge compared to o = 0.98. For the size that we consider here, the agreement is very 
good with the cr c (L) obtained with the domain- wall measurements, see Fig. El Additional 
runs were performed for the value a = 0.98, in particular for larger sizes, up to L = 12. 
For a = 0.98 the average is performed over at least 20 000 samples, each sample being a 
strip of size L x 10 7 . We will employ in the following this data to compute the central 
charge that we will compare to the ±J disorder case. In Fig. ^2 we compare the effective 
central charges obtained for the Gaussian and the ± J distributions of disorder. For the 
±J distribution, we employed the data for p = 0.1095. For comparison, we also show 
the central charge obtained for percolation on an Ising model. In all these cases, we 
show the central charge obtained by a fit to (J4.1)) for increasing L while keeping the 
maximum size fixed to L max = 12. The agreement between the two types of disorder is 
excellent. We see that the central charges converge to ~ 0.465 at the multicritical fixed 
point on the Nishimori line for both cases of disorder. We thus have strong evidence of 
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the universality of the multicritical fixed point. Moreover, the asymptotic value is clearly 
different from the value of percolation, either the numerical one on a finite system or the 
infinite system limit also shown in Fig. ^2 

5 Spin-spin correlation functions 

In the previous sections, we have already seen that the universality class of the multi- 
critical point on the Nishimori line, for the two types of disorder considered, is different 
from the one of percolation. Still the central charge is very close and it is then useful to 
find further measurements to confirm this result. In this section, we present numerical 
results for the spin-spin correlation functions on the Nishimori line close to the Nishimori 
point. 

A general correlation function C n between two points r\ and r 2 has the following 
power-law form with exponent 2 X n 111 £L two-dimensional conformal field theory on the 
infinite plane: 

C n (fi,r 2 )oc— —5^. (5.1) 

\ri-r 2 \ 

In the transfer matrix computations we use long strips with periodic boundary conditions 
along the short directions which one can also interpret as a cylinder. Therefore, we 
consider the correlation function C n on the infinite cylinder of circumference L with 
coordinates u G [1, L] and v e] — oo, +oo[. Applying a conformal mapping to (|5.1j) . one 
infers the following behaviour on the cylinder (see, e.g., chapter 11.2 of [42]) 

C n ((u 1 ,V 1 ),(u 2 ,V 2 )) OC — — — L —— WMI*" • ( 5 - 2 ) 

{2cosh(27r(i>i — v 2 )/L) — 2 cos(27r(wi —u 2 )/L)\ 
There are two particular cases that we will consider in the following: 

• v\ — v 2 , u — U\ — u 2 . The correlation function is measured for two points across 
the strip, separated by a distance u. The correlation function takes the form: 

C n (0,u) oc (sin( — )Lj . (5.3) 

One can then extract the exponent 2 x n by a direct fit of the measured correlation 
function, after averaging over the disorder, to this form. One disadvantage of this 
method is that we can use only points separated by distances u < L/2. 

• When one chooses u± = u 2 , v = v± — v 2 the correlation function takes the following 
form: 

C n (0,v)<x (sinh(^)L)~ 23:n . (5.4) 

Thus we consider two points separated along the strip with the same position across 
the strip, and it is possible to access long distances. 
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We will apply these two forms to the nth moments of the spin-spin correlation function 
[{S(ui, Vi)S(u 2 , V2)) n ], where [• • •] stands for the average over disorder. The general result 
is then, for the two cases discussed above: 

[(S(0)S(x)) n ] oc (sin(^)Lp" , 

[(S(0)S(y)r\ oc (sinh(^)Lp . (5.5) 

Here, we have identified rj n = 2x n . All along the Nishimori line, the moments of these 
correlation functions are equal two by two [2], thus we have the general result that 
Vzk-i — V2k- F° r a pure system, one has i] n = n x rj. On the other hand, in the case 
of percolation over Ising clusters, it is easy to see that all moments of the spin-spin 
correlation functions are equal (and not only two by two). Then, if the Nishimori point 
would be in the percolation universality class, the exponents for the correlation functions 
should collapse to a unique value rj n = r\ at the critical point. 

As a first measurement for ±J disorder, we show in Fig. 1121 the effective magnetic 
exponent rji that we obtain from a fit to the form eq. (|5.3jl . The effective magnetic 
exponent r\\ was computed for three points on the Nishimori line: at p = 0.105, which 
is well inside the ferromagnetic phase and indeed, we see that the effective magnetic 
exponent tends to a small value; at p = 0.115 which is well in the paramagnetic phase, 
which is confirmed by the fact that the effective magnetic exponent increases; at p — 
0.1095 ~ p c where the effective magnetic exponent seems to converge to a value rji ~ 
0.185. On this figure, we also show the corresponding magnetic exponent obtained from 
simulations of the pure Ising model and also of the percolation model, as well as the 
expected values for these two models in the infinite limit (tj — j and ^, respectively 
- see, e.g., [43] and [37], respectively). These two measurements are presented in order 
to show what type of correction we can expect in such a measurement of the magnetic 
exponent. Indeed, we can see that this method gives rather accurate measurements for 
large L (for L = 12, the deviation is around 1%). From this, we can conclude that 
the value of rji for p = 0.1095 is significantly distinct from the one of percolation. This 
is an additional proof that the multicritical point on the Nishimori line is not in the 
universality class of percolation. 

Next, in Fig. we show the same quantity but only close to the critical point. In 
this figure, one can clearly distinguish a change of behaviour close to p ~ 0.1095. For 
lower p, the exponent decreases, which is expected since we are in the ferromagnetic 
phase. On the contrary, for p larger than 0.1095, we observe that the magnetic exponent 
increases, as expected since we are in the paramagnetic phase. 

Fig. El shows the moments of the spin-spin correlation function for L = 20 and 
p = 0.1095. There is a similar plot in [15], but the present data set is completely 
independent. The present data was obtained for a geometry which differs from [15], 
namely on 4001 x 20 strips with a globally fixed number of positive (negative) bonds 
while for that of [15], each bond was assigned a value independently. Here we have 
discarded the 1000 initial (and final) iterations before taking 101 measurements of the 
correlation functions every 20 iterations. The data in Fig. ITU was obtained by averaging 
over 7623 such strips, resulting in statistical error bars which are much smaller than the 
size of the symbols. Despite these differences, the present results for [(S '(0) S \x)) n ] agree 
with those of [15] within error bars. 
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Figure 12: Effective magnetic exponent rji for the ±J RBIM with p = 0.105, p = 0.1095, 
p = 0.115, the percolation model and the pure Ising model. 
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Figure 13: Effective magnetic exponent rji for p = 0.108, • • • , 0.111, for the ± J disorder 
on the Nishimori line. 
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Figure 14: Moments of the spin-spin correlation function for the ±J RBIM with p = 
0.1095 and L = 20. We only show the odd moments. Error bars are much smaller than 
the size of the symbols. 

A direct fit on a doubly logarithmic scale of the correlation functions for p c = 0.1095 
and L = 20 to the form (|5.3|) yields: 

Vl = m = 0.1848 ±0.0003, 

m = r/ 4 = 0.2552 ± 0.0009 , 

V5 = r/6 = 0.3004 ±0.0013, 

r} 7 = r/ 8 = 0.3341 ± 0.0016 . (5.6) 

These estimates are consistent with those of [15,44] within error bars. 

We now turn to the Gaussian distribution of disorder. In Fig. ^1 we present the 
exponent r)i obtained from a direct fit with eq. (J5.3j) in function of a. For that case, we 
have data for size up to L = 14 and for each value of L and a we average over 10 000 
samples of geometry L x 200 L. We clearly see a crossing of the curves close to a = 0.98, 
in agreement with the previous results. At the crossing point, we have rj% ~ 0.18, which 
is very close to the corresponding value for the ± J distribution of disorder. 

Next, we present the exponents obtained, again for L = 14, from a direct fit on a 
doubly logarithmic scale to the correlation functions for a = 0.98: 

Vl = r/ 2 = 0.1818 ±0.0002, 
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Figure 15: Effective magnetic exponent 771 vs. L for the Gaussian distribution of disorder. 

V3 = j] A = 0.2559 ± 0.0002 , 
m = V6 = 0.3041 ± 0.0002 , 

rj 7 = V8 = 0.3402 ± 0.0002 . (5.7) 

The value of these exponents is close to the ones for the ±J distribution of disorder. 
Note that there is still a small difference, r/f J ~ 0.1848 compared to ^Gaussian ^ 0.1818. 
This difference is due to the fact that the exponents are obtained only close to the critical 
points. For the ±J disorder case, the measurement is done at p = 0.1095. From Fig. 1131 
one can read off that this will imply a change of order 0.002 on rjx if we take p c = 0.1093. 
Taking into account this correction, the correspondence of rjx is nearly perfect between 
the two types of disorder, thus giving more support for the universality. 

In Fig. we show the exponent 771 obtained by an independent set of simulations 
on a geometry L x 100 L and with the measurement of the correlation function along the 
strip. We simulated 10 000 samples for each size up to L = 12 with this geometry. The 
exponent 771 is obtained from a fit with eq. (|5.4|) . To perform the fit, we keep only the 
data for correlation functions with two points at a distance y such that 10 < y < 10 x L. 
Thus one does not use the data for two operators very close, contrary to what is done 
while fitting with eq. (|5.3|) and we expect to reduce the finite-size corrections. In Fig. 
one sees that for the largest L, one obtains a constant exponent for a = 0.98 — 0.981 
with rji = 0.180 — 0.182. Thus these results are in perfect agreement with the previous 
measurements on a different geometry. 
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Figure 16: Effective magnetic exponent 771 vs. L for the Gaussian distribution of disorder 
with a = 0.976, 0.978, 0.979, 0.980, 0.981, 0.982, 0.984 (from bottom to top). 
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Figure 17: Thermal exponent vs. L for the ±J distribution of disorder. 

Finally, in Fig. El we present results for the correlation function of the energy opera- 
tor tij = JijSiSj. A thermal exponent v T is associated to this operator, which corresponds 
to a perturbation in temperature, via the relation 



This exponent replaces the exponent x n in eq. (|5.3|) if we replace the spin operator with 
the energy operator. In Fig. El v t is shown for p ~ p c and for increasing strip widths 
L. Finite-size corrections are very strong for ut- We obtain a value of Pt — 3 for 
L — 14, . . . , 16, but it is not clear if we have reached a large enough size. With a different 
method, Merz and Chalker [19] obtained a value of — 4 ± 0.5 by measurements on 
system sizes up to L = 32. We believe that the extrapolation of our results for larger 
sizes is compatible with this result. 

5.1 Central charge and magnetic exponents 

We now return to the central charge. We mentioned at the beginning of section 0] that 
the measured central charge is an effective quantity which can be affected by operators 
with negative dimensions. Indeed, if one considers a torus, i. e. a strip of geometry M x L 
with periodic boundary conditions in both directions, one has the relation [45] 




(5.9) 



26 



with 



Q 



(5.10) 



Here, the index a is associated to the primary operators <ft a which appear in the transfer 
matrix, while the index n& is non zero (and positive) for the conformal descendants only. 
N at b counts the multiplicity of the descendants of the operator <ft a at the level b. Next, 
we consider the free energy per spin: 
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(5.11) 



If all the operators which appear in the transfer matrix have a non-negative dimension, 
then in the limit M 3> L, the last term in eq. ()5.11j) can be dropped and one recovers 
eq. (|4.1|) . On the contrary, if an operator with a negative dimension is present in the 
transfer matrix, then for large M, the last term in eq. (|5.1H1 will be dominated by this 
operator. Let us call 5 < the dimension associated to the operator with the lowest 
dimension. Then eq. (I5.11|) will become, in the limit M ^> L: 
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f(L,M) = f hulk + — (c-12 5) + 



(5.12) 



Thus, if a negative dimension operator appears in the transfer matrix, the measured 
central charge will be only an effective quantity c e s = c — 12 5 > c. If no negative 
dimension operators appear in the transfer matrix, then one has 5 = since the identity 
operator with zero dimension is always present. 

It is known that negative exponents arise at the Nishimori point in certain moments of 
the correlation functions of the disorder operator [20] . Still, these operators will not cause 
any problems in the present context since the disorder operator is not a local operator 
of our theory. 

We know of no direct way to determine if an operator with a negative dimension 
is present in the transfer matrix. Still, one can make the following simple test. In the 
pure Ising model, the lowest dimension corresponds to the magnetic operator. Since one 
expects that the magnetic exponent is present and we know that its dimension is small, 
one can compare it with the first dimension which appears in eq. (J5.11|) . Assuming that 
there is no negative dimension operator in the transfer matrix, the last part of eq. (J5.11|) 
takes the form 



In 1+e 



(-2-4^1) 



ML 
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(5.13) 



(5.14) 



one can estimate directly 5\. In Fig. EE we show this quantity for L = 4, 6 and 8, as well 
as a fit to the data in the range 1.5 < M/L < 5. These bounds are selected by imposing 
a good quality of the fit. The values extracted for 6% are 0.108(1), 0.104(1), 0.102(1) for 
L = 4, 6, 8, respectively. These values seem to converge to a value close to the one of 
the magnetic exponent 771/2 ~ 0.0925 such that it is reasonable to identify it with the 
magnetic operator. This suggests that the measured central charge is indeed the real 
central charge. 
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Figure 18: Fit of the free energy for the ±J RBIM with p = 0.1094 to the form (|5.14j) 
in the range 1.5 < M/L < 5 for L = 4, 6 and 8. 



28 



6 Other measurements 



In this section, we present measurements of other quantities, namely the Binder cumulant 
and the magnetic susceptibility. For each of these quantities, we do not expect to improve 
precision, the measurements are rather done in order to check the consistency of the 
previous measurements. 

6.1 Binder cumulant 

We first present measurements of the magnetic Binder cumulant on the Nishimori line, 
in order to perform an independent measurement of the critical p for the ± J RBIM. 
The magnetic Binder cumulant is defined as follows in terms of the moments of the 
magnetization m [46,47]: 



We note that on the Nishimori line the magnetic Binder cumulant is identical with the 
overlap Binder cumulant which is usually employed in measurements for spin-glass mod- 
els, since we have an equality of the first and second moment of the spin-spin correlation 
functions [2]. 

Simulations are performed on square lattices with periodic boundary conditions in 
both directions and with linear sizes in the range L = 3 up to L = 8. We employ the 
transfer matrix to compute the partition function without and with a small magnetic field 
h as well as with 2h. Then we extract the second and fourth moment of the magnetization 
m from the expansion 



and a similar expansion of Z(2h). The terms (m) and (m 3 ) do not appear in the expansion 
since they vanish at h = on a finite system. We used a value of h ~ 0.01/L 2 in our 
simulations. In order to reach a good precision, a large number of samples had to be 
simulated, typically one million samples for each size and value of disorder p. Fig. ITDl 
shows a plot of the Binder cumulant versus p. In this figure, we observe a crossing in the 
expected region, i.e. p ~ 0.11. Since the number of samples that we have to simulate is 
huge, it is difficult to reach sizes large enough to improve the previous estimate of the 
critical point on the Nishimori line. Thus the measurements that we show here should 
be regarded as a consistency test only. 

The inset of Fig. 1191 shows a plot of the Binder cumulant versus the rescaled variable, 
(p—Pc) L 1 ^. Assuming the values that we obtained in sectional p c = 0.1094 and v = 3/2, 
we see that we have a good scaling behaviour of the Binder cumulant already for small 
lattices. 

6.2 Susceptibility 

A second quantity of interest is the susceptibility \i which we measure using the transfer 
matrix. We compute the free energy, first without any applied magnetic, and next with 




(6.1) 




1 + y (m 2 ) + ^ (m 4 ) + • • ■ 



(6.2) 
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Figure 19: B(L) vs. p on the Nishimori line for the ± J disorder case. The inset shows 
B(L) vs. (p - Pc)L^ u , with p c = 0.1094 and v = 3/2. 
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Figure 20: Susceptibility exponent 7/z/ vs. p on the Nishimori line for the ±J disorder 
case. The inset shows x(L)L~ J / u vs. (p — p^L 1 ^, with p c = 0.1094 and v = 3/2. 
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a small magnetic field. /z,(0) and fh{h) are determined with the same realization of 
disorder in order to decrease the fluctuations. The susceptibility is obtained by the 
following relation 

f L (h) - f L (0) ~ (Ph) 2 X (L) . (6.3) 

We need to choose a very small magnetic field h, typically h ~ 0.000001, in order to 
ensure that we do not have a magnetic term. 

We expect that the susceptibility is of the form x(L) — L 1 ^ at a critical point. In 
Fig. EDI we present the effective susceptibility exponent obtained by two-point fits for 
data with odd sizes. For each of these points, we made 10 000 measurements on strips 
of sizes L x 10 5 . We see that the crossing converges towards a value of p close to 0.109. 
Moreover, we can see that the value of the susceptibility exponent seems to converge 
to a value close to 7/1/ ~ 1.82. Thus one can also check that the hyperscaling relation 
7/1/ + 2/5 /v = dis satisfied with 2(3/v = % ~ 0.18. 

Finally, we can also use the measured values of the susceptibility to perform a fit in 
the following form 

X (L)~L^a((p-p c )LV») . (6.4) 

Since this fit involves three parameters, we obtain a large error bar on each of them. 
Keeping only the reasonable fits, one obtains good collapses of the data in the following 
range: j/v = 1.8 - 1.82, p c = 0.109 - 0.110, l/v = 0.65 - 0.7. In the inset of Fig. EH 
we present a plot of x(L)L~ 7 ^ vs. the rescaled variable (p — p c )L l l u with the values 
p c = 0.1094, l/v = 2/3 and 7/1/ = 1.82 with an excellent collapse of the data. 

7 Dilution 

In this section, we will consider the more general case of a binary distribution with 
dilution. We denote by q the amount of dilution and by p the amount of disorder (see 
section 2 for the definitions). The Nishimori line is now replaced by a surface in the 
T — p — q space. For p = 0, we expect only two critical points [36]: an attractive point for 
q = (no dilution) and a repulsive point for q = 0.5 which is a percolation fixed point 
(see figure EI) and is on the Nishimori surface. Another fixed point, also on the Nishimori 
surface, is the fixed point on the Nishimori line determined previously for q = 0. In 
this section we want to study the flow between these two fixed points on the Nishimori 
surface. In particular, obtaining a clear flow between these two fixed points will give 
further evidence that the multicritical point on the Nishimori line (i.e. without dilution) 
is not in the same universality class as percolation. Our measurements are carried out as 
follows: for a fixed dilution q, we perform simulations for varying p and T in the Nishimori 
surface defined in eq. ()2.8|) and look for a maximum in the effective central charge. Next 
we check the value of this maximum c e g(q) vs. q which is shown in Fig. |^ for central 
charges obtained numerically with a fit to the form eq. (|4.1|) for the sizes [4,5,6,7,8], 
[5,6,7,8] and [6,7,8]. As expected, c e g(q) varies monotonically between c e s(q = 0) and 
Ceg(q = 0.5) showing that there is no additional fixed point. 

In practice, since c e s changes on a very small range (between 0.47769 for q = 0.5 and 
~ 0.464 for q — 0) with strong finite-size corrections, it is more convenient to consider 
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Figure 21: Central charge vs. dilution q along the intersection of the p — q — T Nishimori 
surface with the paramagnetic-ferromagnetic critical surface for the ±J disorder case. 
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Figure 22: X vs. dilution q along the intersection of the p — q — T Nishimori surface 
with the paramagnetic-ferromagnetic critical surface for the ±J disorder case. 
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Figure 23: Magnetic exponent 77 vs. dilution q along the intersection of the p — q — 
T Nishimori surface with the paramagnetic- ferromagnetic critical surface for the ±J 
disorder case. 



the following quantity: 



Ccff - Ceff(0) 



X= 7"^"' (7.1) 
c eff (0.5)-c efr (0) ' 1 ; 

where for each range of sizes L we consider the values obtained with the same corrections 
for c e ff(0) and c e ff(0.5). By construction, this quantity is equal to zero at q = (no 
dilution) and is equal to 1 at q — 0.5 (the percolation case). Fig. |2U shows the quantity 
X as a function of dilution q, again for the sizes [4, 5, 6, 7, 8], [5, 6, 7, 8] and [6, 7, 8]. Apart 
from the size of the errors, one does not observe significant differences. In all cases, one 
immediately sees that the central charge is dominated by the behaviour of the fixed point 
without dilution. If one starts at the percolation point, a small decrease in q will appear 
as a jump towards the fixed point without dilution (q = 0) 5 . At q ~ 0.4, X is already 
indistinguishable from zero, indicating that the central charge is the same as the one for 
q = 0. One expects that by increasing the sizes of the data in the fit for determining c e fr, 
one should observe a crossover between the two fixed points q = 0.5 and q = (such an 
effect will be shown in the next section for the flow between the Nishimori point and the 
fixed point of the pure Ising model). Here we cannot really observe this effect since the 
crossover is too fast. 



5 More precisely, critical properties at finite temperature are expected to be controlled by the fixed 
point N' in Fig. El [4]. Fortunately, the crossover scale for the flow between N and N' appears to be so 
small that this difference can be neglected in practice. 
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Finally, Fig. |2H1 shows the magnetic exponent vs. the dilution measured on long strips 
of width L = 8. Here again, we observe the same effect. For q = 0.5, the magnetic 
exponent is known exactly, it is rj = 5/24 ~ 0.208333 [37]. As we decrease q, the 
magnetic exponent jumps to a value compatible with the one of the multicritical point 
on the Nishimori line 7] ~ 0.18. 

8 Out of the Nishimori line 

In this section, we present some results off the Nishimori line. We will consider two cases 
separately. First we investigate the line which corresponds to the flow from the Nishimori 
point to the pure Ising model. This line corresponds to the ferromagnetic-paramagnetic 
transition line. We show in the next subsection that this line can be determined from the 
maximum of the central charge when varying temperature at fixed p. We show further 
that the flow on this line is from the Nishimori point towards the pure Ising model fixed 
point. Next we turn to the line which connects the Nishimori point with a fixed point 
at zero temperature. An important issue concerns the verticality or re-entrance of this 
line. By arguing that the nature of the Ferro-Para transition at the Nishimori point 
is of geometric origin, it was suggested [3,48] that the transition should take place at 
the same concentration of impurities for any temperature below the Nishimori point, 
implying verticality of the line. However, recent numerical results for the ±J RBIM 
seem to advocate instead a re-entrance of the paramagnetic phase (see Table Q). The 
other important issue is the nature of the fixed point at zero temperature. It was argued 
in light of previous numerical results that the universality class of this point could be 
percolation [17,25], although there is no obvious reasoning supporting this conclusion in 
contrast to the case of dilution. 

8.1 Flow from the Nishimori point to the pure Ising model 

First we study the line which separates the ferromagnetic phase from the paramagnetic 
phase between the Nishimori point and the pure Ising model. The phase boundary itself 
has already been determined with good accuracy [19]. Here we will use a computation 
of the central charge to follow this line in a similar way to what was done in the case of 
dilution. The motivation for studying this line is first to clearly see in which direction 
we flow (we expect towards the pure Ising model fixed point since this is a marginally 
attractive fixed point [1]). Furthermore, we want to check that there is no additional 
fixed point. In a similar study for the 3-state Potts model [49], an additional fixed point, 
predicted by perturbation theory [50,51] was observed. 

We start from the pure Ising model, with a small perturbation, say p = 0.01 and 
T ~ T c , T c being the critical temperature of the pure Ising model. Next we vary T 
and measure the central charge. For small p, these measurements are very simple to 
perform since we have only a weak disorder and moreover, measurements are performed 
at constant p. It is then easy to determine a maximum since the same configurations 
of disorder are employed for different T. Then we iterate the process for larger p and 
follow the ferro-para line which is identified with the maximum of the central charge. 
The measurement is more complicated close to the Nishimori line because the transition 
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Figure 24: T — p phase diagram for the J = ±1 disorder case. 

line has a strong curvature. One needs to change simultaneously p and T, thus we do 
not have any more correlated samples and it is much more time consuming to locate the 
maximum of the central charge (as is the case on the Nishimori line) . Fig. shows the 
phase diagram obtained from the maximum of the central charge. Within error bars, 
our results agree with the phase diagram Fig. 7 determined in Ref. [19] by a different 
method. In Fig. we plot the corresponding central charge obtained from three-point 
fits with L = 3,4,5 and L = 4,5,6. Since the difference of the central charge between 
the pure Ising model and the Nishimori point is very small, we employed a parameter X 
defined as follows: 

x _ C effO) ~ C eS (P = 0-109) (gl) 

c cS (p = 0) - c cS {p = 0.109) ' 1 ' ) 

similar to the one defined in the study of dilution, see eq. (|7.1|) . From this plot, we can 
see that by increasing the lattice size, the attractive fixed point is the pure Ising model 
(since the variation of the central charge increases, starting from the Nishimori point). 
This is consistent with a flow from the Nishimori point to the pure Ising model and 
moreover the absence of additional fixed points along this line. 



8.2 The zero-temperature fixed point 

In this subsection, we study the zero-temperature fixed point. First we present results for 
the location of the fixed point at T = 0. This was already determined with high precision 
for the ±J disorder case by Wang et al. [12] and Amoruso and Hartmann [28]. The result 
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Figure 25: Variation of the effective central charge as characterized by the quantity X 
(see (J8.ip ) on the Ferro-Para line, for the J = ±1 disorder case. 
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obtained by these groups is p° c ~ 0.103, indicating a re-entrance of the paramagnetic 
phase. For comparison, we repeated these measurements, obtaining compatible results. 
We further considered the Gaussian case for which the same re-entrance effect is obtained. 

The simulations were performed on square lattices with free boundary conditions 
along one direction and periodic/antiperiodic conditions in the other direction. We mea- 
sure the difference of energy between the periodic case and the antiperiodic case 

AE = E p — E a . (8.2) 

This difference of energy corresponds to the energy associated to a domain wall in the 
system. We expect that after averaging over disorder, the domain-wall energy is charac- 
terized by some size exponent p 

[AE] oc L p . (8.3) 

This quantity (as well as the one associated to the width of the distribution of domain- 
wall energies) was computed with high precision by Amoruso and Hartmann [28] on very 
large lattices, up to L = 700, by using a minimum-weight perfect matching algorithm. 
In Fig. I2H1 we present our results on smaller sizes (up to L — 100) but with much larger 
statistics, while employing the same type of algorithm [52]. We simulated 1 million 
samples for each p and L < 100 and 0.5 million samples for L = 100, compared to 30 000 
samples in [28]. The reason for desiring even better precision is to be able to observe a 
crossover in the location of the fixed point, see the discussion of this point in the following 
subsection. In Fig. |2H1 we clearly distinguish a scaling in function of the size L close to 
p ~ 0.103 for large L, in perfect agreement with the previous results [12,28]. We also note 
that this power- law behaviour is apparent only for sizes L > 20. In Fig. we present 
the same quantity for the Gaussian distribution of disorder. Here again we observe a 
power law of [AE] with size L close to er ~ 0.97 to be compared to the result (|3.11|) 
a c ~ 0.97945 on the Nishimori line. For the Gaussian disorder, we have data up to size 
L = 140 and the number of samples is again 1 million for each value of a and L except 
for L = 140 where we have 0.5 million samples. 

Thus in both cases, we observe a re-entrance of the paramagnetic phase. 
8.3 Magnetic exponent at zero temperature 

In order to characterize the fixed point at zero temperature we measure the magnetic 
exponent in the same manner as on the Nishimori line. Specifically, we use measurements 
of the spin-spin correlation functions which provide a direct estimate of the magnetic 
exponents and allow us to study higher moments. 

Fig. EH] shows the values of rji, obtained from the measured spin-spin correlation 
functions with a fit to the form eq. ([5.3)1 . We see in this figure that rji is constant close 
to p = 0.107 — 0.108, at least for the largest sizes that we can reach, L = 20. For larger 
sizes, we expect that this value can still decrease, see the discussion below. Next, in Fig. 
|2*9"1 we compare the magnetic exponent (or more precisely rji = 2xh) obtained for the ± J 
distribution of disorder, both on the Nishimori line and at T = 0. In this figure, one can 
see that the critical point at T = 0, denoted by p° c is very close to p c , the critical point 
on the Nishimori line. The best measurements yields a value of p° c slightly smaller than 
p c , close to 0.108. This value is far from the one obtained in the previous section, i.e. 
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Figure 26: Domain-wall energy [AE] as a function of the system size L for the ±J 
disorder at T = and for p = 0.100, . . . , 0.105. 
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Figure 27: Domain-wall energy [AE] as a function of the system size L for the Gaussian 
disorder at T = and for a = 0.96, . . . , 0.98. 
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Figure 28: Effective magnetic exponent r]\ vs. L for the ±J disorder case at T = 0. 

p° c = 0.103. We explain this difference by the existence of strong finite-size corrections, 
a situation which is frequent for two-dimensional spin glasses, see [53]. These finite- 
size corrections can also be observed directly in Fig. |2H For small sizes, 6 < L < 10, 
the effective exponent rji is constant for p ~ 0.109 and only for this value of p. For 
10 < L < 14, it is almost constant for p ~ 0.108. By increasing L, the estimate for p° c 
will continue to decrease. For the largest size that we simulated, p° c is estimated close to 
0.107. Presumably, this value will still decrease with increasing size. In the discussion in 
the previous subsection, we had already observed that a power-law behaviour is apparent 
only for L > 20. 

Fig. EDI shows 7/1 for the Gaussian distribution of disorder on the Nishimori line and 
at T = 0. Here again, we clearly observe a small re-entrance, with cr° ~ 0.97 compared 
to a c ~ 0.97945 (see ()3.11|) ) on the Nishimori line. Contrary to the ±J disorder case, 
there are very small finite-size effects. For L ~ 14, one already obtains the same result 
for (7° as by domain-wall measurements on much bigger systems. 

Finally, let us compare the value of the magnetic exponents for the two types of 
disorder at T = 0. For the Gaussian disorder, one obtains r/i = 0.11 — 0.12, see Fig. EH 
The ± J disorder case is less clear. In Fig. |2H as discussed above, one still has strong 
finite-size corrections. If one considers the asymptotic point p\ ~ 0.103 (compare section 
18.2)) . then the corresponding rji will be close to 0.11, thus rather close to the one of the ± J 
disorder case. Still, this is not sufficient to conclude that we have universality at T = 0. 
Indeed, at T = 0, the type of disorder influences the degeneracy of the ground states. 
For the Gaussian distribution of disorder, each correlation function (before averaging over 
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Figure 30: Effective magnetic exponent T]i vs. cr for Gaussian disorder on the Nishimori 
line and at T = 0. 
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disorder) will be equal to ±1. Thus the even moments will always be equal to one, while 
the odd moments will be equal to the first moment. Then, one obtains the general result, 
valid as long the ground state is unique (up to a global symmetry), r] 2n+ i = T]i for any n 
and rj2n = 0. On the contrary, for the ±J disorder case, we have a huge degeneracy of 
the ground states. A direct determination of the exponents rji for L = 20 and p = 0.107 
yields: 

77! = 0.1277 ±0.0004 ; rj 2 = 0.0593 ± 0.0002 
773 = 0.1369 ±0.0004 ; 774 = 0.0719 ± 0.0002 
775 = 0.1412 ±0.0004 ; 77 6 = 0.0776 ± 0.0002 
777 = 0.1438 ±0.0004 ; % = 0.0810 ± 0.0002 . (8.4) 

The value p = 0.107 is chosen even if it is not exactly at the asymptotic critical point 
p° c ~ 0.103 since for this value we have measurements for large sizes. We also computed 
the same exponents for smaller p and smaller L. The change in the exponents is very 
small. We thus conclude that for the ±J disorder case, all moments are different. 

9 Summary and conclusions 

In this paper we have performed an extensive study of the random-bond Ising model 
with a particular emphasis on non-trivial fixed points. Concerning the Nishimori point, 
which is located at the intersection of the Ferro-Para critical line and the Nishimori line, 
we have used domain-wall free energy computations to provide an accurate estimate for 
the location of this multicritical point. This has been done for both, the binary ±J 
distribution as well as for the Gaussian distribution. Our numerical results show that a 
conjecture for the location of the Nishimori point based on a duality argument [21] yields 
only a very good approximation, but is not exact. 

Next, we estimated v = 1.48(3) for ± J disorder and v = 1.52(3) for Gaussian disorder. 
This agrees with other recent estimates for ± J disorder, namely v = 1.50(3) on the square 
lattice [19], and v = 1.49(2) on the triangular and honeycomb lattices [54]. All these 
results are consistent with a universal value v « 1.50 for the Nishimori point. We have 
obtained further accurate results for the exponents 771, . . . , rjs for the moments of the spin- 
spin correlation functions. The estimates for ±J disorder (|5.6|) and Gaussian disorder 
(15. 7|) are not only very close to each other, but also to recent estimates on the triangular 
and honeycomb lattices [54] . Our analysis of the central charge c at the Nishimori point 
for both types of disorder is also consistent with a universal value c = 0.464 ± 0.004 [15]. 
All these results 6 suggest a single universality class of the Nishimori point in the two- 
dimensional random-bond Ising model, and definitely exclude percolation as the possible 
universality class for this point. 

We have also considered a probability distribution for the bonds with dilution, in 
which some of the coupling constants are zero. In the purely diluted case (a distribu- 
tion containing only positive or zero bonds) there is only one non-trivial fixed point, 

6 Essentially the same results are found [55] if the approximate value p c m 0.110028 [21] is used instead 
of the numerically exact value l|3.4|l . 
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namely the zero-temperature percolation point, apart from the critical point of the pure 
system [36] . We have confirmed that this percolation point is unstable against the Nishi- 
mori point if one moves within the intersection of the critical transition surface and the 
Nishimori manifold by considering both, dilution and ± J couplings, and that these two 
points (percolation and Nishimori) are the only fixed points within this intersection line. 
On the other hand, going off the Nishimori line but staying within the critical transition 
line we confirm that the Nishimori point is unstable in favour the pure Ising model fixed 
point. All these results are obtained by studying the crossover of the effective central 
charge, or first Lyapunov exponent, in strips of increasing width. 

Finally, we have analyzed the critical point at zero temperature corresponding to the 
ferro-para transition in the model without dilution for both the binary and Gaussian 
distributions. Our numerical analysis confirms the strict re-entrance of the ferromag- 
netic phase [12,28]. We have also investigated the criticality of the zero-temperature 
fixed point and argued that it is, as for the finite-temperature Nishimori point, different 
from percolation. The results obtained in this paper raise the question of the apparent 
hierarchy of fixed points for strongly disordered systems, namely in our case the pure 
Ising transition, percolation, the Nishimori point and the zero-temperature point, all cor- 
responding to different (and certainly non-unitary) conformal field theories with values 
for the central charge and critical exponents which are extremely close to each other but 
nevertheless different. Future analytical efforts will be needed to understand such effects 
of strong disorder in two-dimensional classical statistical systems. 
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